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O ' Abstract 
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New type of tomographic probability distribution which contains 
complete information on the density matrix (wave function) related 
to the Fresnel transform of the complex wave function is introduced. 
O ■ Relation to symplectic tomographic probability distribution is eluci- 

dated. Multimode generalization of the Fresnel tomography is pre- 

Ch ' sented. Examples of applications of the present approach are given. 

d ' 
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1 Introduction 

Recently symplectic tomography introduced in quantum òptics [T] and sig- 
nal analysis j3 EI was extended to the problem of soliton solutions for some 
nonlinear dynamical equations like nonlinear Schròdinger equation jH |3] in- 
cluding the study of solitons in Bose-Einstein condensate jü] described by 
Gross-Pitaevsky equation jjj. The name "symplectic tomography" is re- 
lated to using in the approach Standard symplectic group transformations in 
phase space (linear canonical transformation in position and momentum pre- 
serving Poisson brackets in classical case and commutator in quantum case). 
The symplectic tomographic map associates with a function (analytic signal, 
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soliton solution, etc.) a specific probability distribution function (symplec- 
tic tomogram) containing the same information on the signal as the initial 
function. It was found in |U |ü] that due to the structure of symplectic 
tomogram depending on three real coordinates, one can construct another 
probability distribution (Fresnel tomogram). This tomogram depends on two 
real variables but information contained in the Fresnel tomogram is sumcient 
to reconstruct the initial function for which this tomogram is known. The 
inverse formula which permits to make explicit reconstruction of the initial 
function (signal, soliton, etc.) was not written until now and it is one of the 
aims of this work, namely, to present the reconstruction formula. Another 
goal of the work is to extend the Fresnel tomography to two-particle (mul- 
tiparticle) systems to be able to consider in the tomographic representation, 
e.g., soliton solutions of nonlinear equations for 2D and 3D cases. 
The paper is organized as follows. 

In Section 2, we review the symplectic tomography scheme both for 1D 
case and multidimensional case. In Section 3, we present the bàsic concept 
of the Fresnel tomography for 1D case while in Section 4 we present an 
extension of the Fresnel tomographic approach to multidimensional problems. 
In Section 5, we apply our method to the case of 1D chirped Gaussian wave 
functions. Finally, the conclusions are summarized in Section 6. 



2 Symplectic tomograms 

Below we present symplectic tomograms associated with a complex func- 
tion ip(x) describing either a quantum normalized state or analytic signal 
or soliton solution of a nonlinear equation. The tomogram is given by the 
formula |2, 



w(X, /i, v) 



2ttM 



^(í/)exp 



iX 



v~ y) d v 

2v v 



(1) 



One can see that the tomogram which depends on three real variables has 
the homogeneity property 



w(XX, À/i, \u) 



|A| 



w(X, //, v). 



(2) 
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The inverse transform which relates the tomogram to the complex function 
(density matrix) reads [S] 



1>{X)r(X') = ±- /«;(Y>,X-X')exp 

Z7T J 



i Y 



I' ■ 



dfidY. (3) 



Since one can reconstruct the density matrix in view of the known sym- 
plectic tomogram, all other characteristics like Wigner- Ville function [ÜJ EI] 
can be also expressed in terms of the symplectic tomogram fTJ 

dX dfi du 



W^(q,p) = w(X,fi,v)exp i(X - fiq - vp) 



(27T) 5 



(4) 



where 



u 



W^(q, p) = — / V ( q + 77 ) r ( Q ~ tj ) e ~^- 



The Wigner function is normalized W^(q,p) dqdp = 1. The symplectic to- 
mogram provides the possibility to find the optical tomogram simply taking 
the parameter of the tomogram in the form /x = cos#, z/ = sinó 1 . One arrives 
at the formula for the optical tomogram in the form 

2 



w(X, 



2tt \ sin 9\ 



if>(y) exp 



i cot 9 



V 



iX 
sin 9 



V dy 



(5) 



Thus optical tomogram can be used in order to reconstruct Wigner func- 
tion fHlE] by means of Radon transform [15 . 

The generalization of symplectic tomogram to the multidimensional case 
is straightforward. 

To provide this generalization we introduce the following notation. 

There are three real vectors 



X — (Xi, X 2 , . . . , X N ), fl — (fií, /i 2 , 



v 



and all the components of these vectors vary from — oo to oo. The tomogram 
w(X, fl, v) associated to a complex function if>(x), x = (xi, x 2 , ■ ■ ■ , x^) reads 

\ 2 

i l· l j „,2 



w(X,fl,v) 



1 



' N 



J j=l \ zz/ j 



iX 3 



y.J 



(6) 



The inverse formula has the form 



1 \ N r - - - í N 

-) / w (y, A x-xo(|Texp 



M Y k - fji k 



X k + XI 



dfïdY. 



(7) 
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3 Fresnel tomography in 1D case 



In jUHniHTj the Fresnel tomogram wp(X, v) was introduced 



w F (X, v) 



V2 



1Ï11S 



exp 



i(x- y y 



2u 



■fp(y) dy 



(8) 



Another version of Fresnel tomography based on using the Fresnel integral 
to reconstruct Wigner function of quantum state was suggested in [TH] . 

Since the symplectic tomogram w(X, //, v) is homogeneous function, there 
exists the relation 



Wf 



X v_ 



(9) 



This means that if one knows Fresnel tomogram wf(X',u'), the symplectic 
tomogram is obtained by the substitution of the arguments of this function 



X' 



X 

n 



v 



v 



and multiplication of the obtained function by the factor The formu- 

lated relations provide the inverse formula to reconstruct density matrix by 
expressing it in terms of Fresnel tomogram 



Hx)r(x') 



i 

2^ 



1 



'Y X-X' s 



x + x ,s 



wp [ — , ) exp i ( Y — fj, ) dfi dY. 

(10) 

One can also reconstruct the Wigner function expressing it in terms of Fresnel 
tomogram 



W{q,p) 



1 



1 (X v 

■ wf 



[i /i 



exp [— i(fiq + vp — X)} dfi dv dX. (11) 



The inverse formula obtained gives a possibility to find the Wigner function 
and the density matrix in the position representation if one knows the Fres- 
nel tomographic probability distribution wp(X, v) > 0, which satisfies the 
normalization condition 



J w F (X, v) dX = l. 



(12) 
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4 Fresnel tomogram for multipartite system 



Let us introduce the Fresnel tomogram for function ip(x) depending on several 
variables. This tomogram can be defined in terms of symplectic tomogram, 
i.e., 

Wp(X, v) — w(X, jl — 1, £7), (13) 

where we have symplectic tomogram and /2 = 1 means that the vector jl is 
taken to have all the components equal to one. Thus the Fresnel tomogram 
is expressed in terms of the complex function ip(x) as follows: 



w F (X, v) 



N 

n 



(14) 



To get the inverse of the above transform, we have to use the multidimen- 
sional analog of relation (jïïj) which reads 



w(X, jl, v) 



Wp 



(15) 



|/il/i 2 • • • fJ,N\ 

where vectors in the argument of Fresnel tomogram have the components 

X\ x 2 x N 

) 5 • • • 5 

Ui /i 2 



Vu 



^1 ' (*2 



Un, 



Using relation (|15|) one can construct the expression of the density matrix in 
terms of the Fresnel tomogram analogous to 1D expression (JTÜj) 



^x)r(x') 



djl dY 
(2^p 



Wf 



N 



x 4 



X': 



(16) 



5 The chirped Gaussian case 

In order to verify the validity of the method, we have numerically simulated 
the reconstruction of a 1D complex Gaussian chirped field (GCF) given by 
the following normalized form: 

/ 2 \ 1/4 

^( X ) = ( ^2 ) eXP 



-— + lax 
a 2 



(17) 
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having a chirp a and a width determined by the parameter a. The above 
state is also known as correlated coherent state [19... According to Eq. (J3J) 
the complex field given by Eq. ()17j) can be retrieved from its tomographic 
representation w(X,fi,x) by the following inversion integral: 

V>O)^*(0) = ^- J J w(X,fx,x)exp[i(X -//x/2)] dX dfi. (18) 

The GCF tomogram can be calculated from Eqs. and ()17j) . The 
Gaussian distribution can be integrated and the tomogram can be written in 
the following form using, according to the notation of the previous sections, 
the symbol v instead of x 



w(X,fM,u) 



2a 2 { 2a 2 X 2 

exp 



\\ 7T [Au 2 + a 2 (fi + 2avf\ "~ r { 4z/ 2 + a 4 (/i + 2avf . 



(19) 

Equation ()19|) shows that the GCF tomogram is still characterized by a 
Gaussian distribution law but its width u is, in general, different from the 
width of the original GCF and it can be expressed in terms of the phase-space 
variables fi and u, chirp parameter a, and width a, namely, 

, _ (/i + 2^)V + 4^ 
U ~ 2^ • (20) 

It can be easily verified that the tomogram w(X,fi, v) in Eq. (fTÜj) coin- 
cides with the tomogram of a Gaussian field without chirp (a = 0) computed 
for the vàlues of the phase-space variables fi = fi + 2az/ and v = v. We 
can admit a different transformation of the phase space, taking into account 
that the tomographic representation of the GCF can be also obtained from 
that of a Gaussian field without chirp with fi = fi and v = vfij^fi + 2av). 
In both instances, the effect of the chirp is that of shifting one of the two 
vàlues of the phase-space variables while leaving unaltered the other one. To 
illustrate the behaviour of the GCF tomogram, we have displayed in Fig. 1 
the 3D distribution of the tomographic map wiXj 'fi,l,v / 'fi) for increasing 
vàlues of the chirp parameter. This representation of the tomogram in terms 
of two parameters X/fi and vj fi rather than three is quite convenient since 
it allows one a clear visualization of the tomographic map in the 2D phase- 
space plane [Xj fi, v / fi). Indeed, recalling the general homogeneity property 
of the tomogram we can obtain complete dependence of the tomogram from 
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Figure 1: 3D plot of the tomogram distribution w(X/fi, 1, v j \i) of the Gaus- 
sian chirped wave function for a width of the Gaussian a = 1 and for increas- 
ing vàlues of the chirp parameter a = 0.5 (a), 1 (b), 2 (c) and 3 (d). 

its three space variables by virtue of Eq. Q : 



According to the Fresnel-based interpretation of the optical tomogram 
discussed in the previous section, w(X//j l , l,u/fj,) represents the GCF inten- 
sity at distance v and Eq. (J21)) telis us that the general dependence of the 
tomogram of its three real variables X, fi, and v can be recovered from a 
set of measurements of the intensity distributions of the propagated GCF 
performed at different distances with a varying scale factor Equa- 
tion (|21|) can be verified immediately from the defining expression of the 
GCF tomogram given by Eq. (|19|). Figure 1 shows clearly that tomographic 
distributions shrink in the phase space with increasing the chirp parame- 
ter of the Gaussian field. The numerical results have been obtained for a 
GCF of width a = 1. In Fig. 2 we display the plots of the tomograms for 
the same set of vàlues of Fig. 1 except that now the width of the GCF is 
a = 0.5. As can be seen by comparing Figs. 1 and 2, the shrinking of the 
3D tomographic distribution in the phase space with increasing the chirp 



W(X, yU, V) 



1 



w(X//i,l,v/n). 



(21) 
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Figure 2: 3D plot of the tomogram distribution w(X/fi, 1, z///i) of the Gaus- 
sian chirped wave function for a width of the Gaussian a = 0.5 and for 
increasing vàlues of the chirp parameter at = 0.5 (a), 1 (b), 2 (c) and 3 (d). 

parameter is less pronounced when reducing the GCF width. A computa- 
tionally efficient method for retrieving the wave function makes direct use 
of the inversion integral given by Eq. (|18p which allows one to determine 
ijj(x) up to the complex constant quantity ^(0)*. The method employs the 
fast Fourier transform (FFT) algorithm for computing the two-dimensional 
Fourier transform w{u x ,u ÍJi , v ) °f the tomogram, namely, 

] /*oo /*oo 

w(w x , uj^v) = — / w(X, /i, v) exp [i (u x X + uj^/j,)} dXdfi, (22) 

Z7T J— oo J— oo 

where u x and u; M are the spatial freqüències corresponding to the phase-space 
variables X and fi of the tomogram. The inversion integral is a particular 
case of the two-dimensional Fourier transform of the tomogram given by Eq. 
(J22|) . In fact, we have 

i/>{u)i/>(Qy = w(l,-v/2,v). (23) 

Equation shows that we can determine the complex quantity ip{y)ip{Q)* 
from the vàlues of two-dimensional Fourier transform of the tomogram at 
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freqüències oox an d lü^ where v is allowed to vary in the range along which 
the tomographic measurements are performed. In principle, this procedure 
works well for reconstructing a wave field, if we have a sufficient number of 
sampled vàlues of the tomographic distribution w(X,n,u), i.e., if the sam- 
pling rate is at least larger than twice the Nyquist rate. In this case, the 
discrete two-dimensional Fourier transform w(rAux, sAu^) corresponding 
to the continuous Fourier integral given by Eq. (J22)) can be written in the 
following form: 

N—l N-l r / x 



- 2n l N + W) 



w(rAüüx, sAujfj) = AJA/i ^ ^ w(nAX,mA/i, i/)exp 

n=0 m=0 

(24) 

where r, s and n, m are integer numbers, AX and A/i are the sampling inter- 
vals in the (X, /x) phase space and Alüx and Au^ are the corresponding sam- 
pling intervals in the frequency space. Equation ()24)1 allows one to compute 
a matrix of iV 2 complex numbers w(r Aux, sAu^) from the sampled vàlues 
w(nAX, mAfi, v) of the tomogram of the wave function via the discrete two- 
dimensional fast Fourier transform (FFT) algorithm. It should be remarked 
that in order to obtain accurate reconstruction of the wave function ^(v), the 
two-dimensional FFT algorithm needs to be repeatedly applied to every two- 
dimensional tomographic distribution corresponding to each v value in the 
considered range and, if high accuracy is required, even this FFT-based re- 
construction method becomes computationally expensive. However, in many 
cases of interest, symmetry consideration allows one to simplify somewhat 
the task through the consideration of a reduced set of sampled vàlues of 
the tomogram. In the GCF case, the two-dimensional Fourier transform 
w{üJx, ^tii v) can be readily obtained in the following form: 



w{u x ,v li ,i') = \ — 2~~r ex P 



x 



v 2 fuj x \ 2 2 / '·· \ 2 

O ) O* \iüx 



2 , ■ - 2iauj fl is 



(25) 



By using Eq. (J2ljj) it can be easily verified that Eq. (|2^|) gives correctly 
■^(^^(O)* with the wave function ip(v) given by the Gaussian chirped field 
in Eq. JT7J). 

6 Conclusions 

We summarize the main results of the paper. 
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We obtained the inverse formula which in explícit form reconstructs the 
wave function, density matrix and Wigner function from the known Fresnel 
tomographic probability distribution function. The Fresnel tomogram can 
be used to measure the phase of radiation by measuring the intensity of the 
radiation. 

We introduced Fresnel tomogram for multimode system and found re- 
construction formula for wave function and density matrix in terms of the 
tomograms for this case too. An example of soliton solution used also in 
analyzing some states of Bose-Einstein condensate can be studied from the 
viewpoint of Fresnel tomography |ïï] . 

The symplectic tomography [T] has as a partial case the optical tomogra- 
phy procedure. In this work, we clarified that there exists another important 
partial case of the symplectic tomography which is named Fresnel tomogra- 
phy. This mean that there exists a connection of Fresnel tomography with 
well-developed optical tomography scheme. The tomographic probability of 
optical tomography scheme is connected with Fresnel tomogram by the rela- 
tion: 

1 í X \ 

w(x, cos 6, sin 6) = — wp ( t , cot 6 ) . 

| cos 1 V sin ) 

One can generalize this formula to multidimensional case. 
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